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Abstract 

We investigate the magnetization reversal processes on classes of complex spin net- 
works with antiferromagnetic interaction along the network links. With slow field 
ramping the hysteresis loop and avalanches of spin flips occur due to topological in- 
homogeneity of the network, even without any disorder of the magnetic interaction 
[B. Tadic et al, Phys. Rev. Lett. 94 (2005) 137204]. Here we study in detail proper- 
ties of the magnetization avalanches, hysteresis curves and density of domain walls 
and show how they can be related to the structural inhomogeneity of the network. 
The probability distribution of the avalanche size, N s (s), displays the power-law 
behavior for small s, i.e. N s (s) oc s~ a . For the scale-free networks, grown with pref- 
erential attachment, a increases with the connectivity parameter M from 1.38 for 
M = 1 (trees) to 1.52 for M = 25. For the exponential networks, a is close to 1.0 
in the whole range of M. 

Key words: antiferromagnets, Barkhausen noise, evolving networks, Ising model, 
Monte Carlo simulations, simple graphs and trees, spin dynamics 
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1 Introduction 



Growing networks is a vivid area of interdisciplinary sciences, with long list of 
applications [1,2,3,4,5,6]. The mainstream of literature is concentrated on the 
structure of the networks. However, from the point of view of some applications 
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it is of interest to decorate nodes with additional degrees of freedom. In the 
simplest case, these variables are discrete, with two possible values: ±1. Such 
an extension can be useful when discussing numerous examples, from neural 
and logical networks through sexual networks to quantum gravitation [7]. Our 
topic here is a network of interacting magnetic moments (spins), with two 
possible magnetic orientations, up or down. The models of spin networks with 
ferromagnetic interaction have been considered in several recent studies [8,9]. 
The case of spin networks with antiferromagnetic interactions is essentially 
different because of the frustration effect along closed cycles on the graph 
[10,11]. This difference is well known also for regular lattices [12]. 

Recently we introduced and investigated field-driven dynamics of spins on 
complex networks with the antiferromagnetic interaction [11]. We found that 
the structural complexity of the networks leads to avalanches of spin flips and 
a criticality of the hysteresis loop. The desired shape of hysteresis curves can 
be obtained by tuning the clustering parameter M. This parameter counts the 
number of links, which attach each new node to the growing network. 

In this work we expand the study of the reversal processes in the antiferro- 
magnetic spin networks. We demonstrate the genesis of avalanches due to in- 
homogeneity in network's connectivity when the clustering is low. In this limit 
the avalanche structure and duration can be clearly interpreted in terms of 
the theoretical distribution of connectivity and depth of the network. For large 
clustering, the character of the avalanche spectra N s (s) for large s change from 
mulitipeaky to a soft decreasing. For small avalanches, the avalanche spectra 
display the power-law character in the whole range of the clustering parameter 
M. We define the size s of avalanche as a number of spins which were flipped 
at a given value of external field H. 

In Section 2 we introduce details of the model and discuss spin reversal on 
small graphs. Section 3 contains details of the emergent structures for two 
classes of growing networks exponential and scale-free. Our results on the 
magnetization reversal in these networks are given in Section 4 and discussed 
in Section 5. 



2 Spin reversal on small graphs 



To describe our procedure in details, we begin with small simple graphs (Fig. 
1). A spin Si = ±1 is placed at each node. If two nodes are linked together, 
antiferromagnetic interaction tends to keep them in one of two antiparallel 
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Fig. 1. Toy graphs and their hysteresis loops for decreasing field (continuous red 
line) and increasing field (dotted black line). 

configurations. Energy of the whole system is written as 

E=-J 1 £S i S j -H y £S u (1) 

(ij) i 

where coupling constant J = — 1 for the antiferromagnetic interaction and the 
first summation goes over all nearest neighbor pairs (ij). We assume symmet- 
ric neighbors, i.e. if spin i is neighbor to spin j then also j is neighbor to i. (We 
note, that results may be drastically different when the neighbor relations are 
directed [9,13].) To draw the hysteresis loops, graphs are placed in saturating 
magnetic field H > 0, which forces all spins to be positive, i.e. to point up- 
wards. We ramp the field H by integer steps [14] from n + 5 to n — 1 + 5 in order 
to avoid ambiguous spin orientation at integer field values H = n (n e Z) . As 
the field decreases gradually at zero temperature T = some spins are flipped 
[Si — > — Si) because of the antiferromagnetic interaction with their neighbors. 

In Fig. 1 we present some graphs and their hysteresis loops. We note that for 
some of them, the obtained hysteresis loop depends on the order of updating. 
As an example, consider a chain of five spins presented in Fig. 1(c). In satu- 
ration, effective fields at node 2, 3 and 4 are the same. If the central spin 3 
is updated at first, only this spin flips at field H = —IJ . If we update spins 
starting from the end of chain (as it is done here), the flipped ones are 2 and 
4. Then, not only the flipped spin labels but also the number of flipped spins 
is influenced by the order of updating. However, any procedure suffers this 
ambiguity until we do not add a small variation of the magnetic moments or 
the interactions; in this case, the Pardavi-Horvath algorithm [15] is appropri- 
ate, where the field is changed continuously and a spin where the effective 
field is zero is flipped as first. On the other hand, the synchronous scheme of 
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Fig. 2. (Color on-line). Snapshots from simulation for scale-free tree (M = 1, 
left) and scale-free simple graph (M = 2, right) for N = 200 and magnetic field 
H = — 5. Spins Si = —1 (red) and Si = +1 (yellow) and their pairs S{Sj = +1 
(black) SiSj = -1 (blue) for all 1 < i,j < N. (Figures using Pajek [21].) 

spin updates may lead to persistent oscillations as when — > —Si, and the 
avalanche will never stop. Here we apply causal way of spin updating, i.e. from 
the spin decorating "the oldest" node of the graph towards the one assigned 
to the node added at very end. Below we present results averaged over a large 
number N run of graphs with different shapes, and the influence of the order of 
updating is expected to be averaged out. 



3 Structure of networks 



It is well known that the structure of the growing network emerges in the 
process of node adding and linking. Here we are interested in the growing ex- 
ponential (EXN) and the scale- free (SFN) graphs [4]. In both network 
grows by successive adding of new nodes. Each new node is linked to M differ- 
ent pre-existing nodes, which are selected randomly (EXN) or with preference 
(SFN). Here we limit our interest to the case when the preferential probability 
of selecting a given node is proportional to the degree k of this node, i.e. to the 
number of currently present bonds at that node. This procedure leads to the 
scale-free character of the degree distribution for SFN, i.e. P(k > M) oc k~ J , 
where theoretically 7 = 3.0 for large network size iV independently on the 
parameter M [3,4,5,16]. For EXN trees (M = 1), P(k > M) = 1' k [4], and 
for EXN with M = 2 we get P(k > M) = 3/4 • (3/2)~ fe [17]. The parameter 
M can be tuned to modify the sparseness of the networks, from M — 1 (trees) 
to higher values, where cyclic paths appear. If the number of nodes in such 
a cycle is odd, as in a triangle, a frustration of spins occurs: some antiferro- 
magnetic bonds join spins of the same orientation, and the minima of energy 
become shallow. The clustering coefficient C [3] depends on the node age. The 
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clustering profiles, as shown in [11] for the case M — 5, exhibits a power- law 
tail both for EXN and SFN. The networks with low clustering (M = 1,2) 
reveal effects of linking inhomogeneity. 

We stress that all bonds are antiferromagnetic and there is no bond disor- 
der in the traditional sense. However, the structure of the network cannot be 
treated as ordered in the sense of periodic lattices. In most of graphs each 
node has its unique position in the network [18]. Therefore, a question arises, 
if the spin glass phase is possible in the antiferromagnetic spin-networks with 
clustering M > 1 at least at T = 0. In terms of the accepted spin-glass termi- 
nology, the presence of nodes with large connectivity may be regarded either 
as a long-range interaction effect, or as an equivalence of a large system di- 
mensionality. However, such structure is not translationally invariant [19]. The 
network shows a nontrivial connectivity profile, which implies new features of 
the spin-glass order and requires refined theoretical approaches [20]. 

In Fig. 2 examples of the networks considered in this work are presented. For 
clarity of the plots the network size is kept moderate. In the simulations we 
use networks with N = 10 4 nodes. 



4 Hysteresis loops and avalanche distributions 

Results of the dynamics of the remagnetization process are shown in Figs. 3 
and 4 for EXN and SFN, M — 1, 2. In Fig. 3, examples of the hysteresis loop 
are presented for EXN and SFN, M — 1,2. The loops are driven in the applied 
field range from —H mSuX — 1 to iJ max + 1, where H mSuX is the saturation field 
equal to the maximal connectivity k max of the lattice. All loops show jumps 
of magnetization at integer field values, as those in Fig. 1. 

In Fig. 4(a) we show the data on the number of flipping spins against the 
magnetic field, calculated within the same numerical experiment, as in Fig. 3. 
In descending field, all the plots in this figure show an asymmetry with respect 
to the zero field. Most of spins flip at H — —1 for trees, H = — 2 for the case 
M — 2. This is a consequence of the connectivity, which is two times smaller 
for trees in the average. 

The results shown in Fig. 4(b) are equivalent to a plot of the interaction energy 
against field. Each antiferromagnetic pair (AFP) gives the contribution J to 
the energy. Then the energy of interaction can be written approximately as 
J (2 N afp ~ MN). Again we see that higher connectivity enhances the amount 
of frustration. We note that the width of the presented curve is much larger 
for SFN than for EXN. This is a consequence of the fact that a node with 
maximal connectivity has much more links for the scale-free graphs. 
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(a) scale-free, M.1, N=10 4 , N run =1 (b) scale-free, M=2, N=10 4 , N m „=1 




Fig. 3. The field dependence of the magnetization m(H). The field is expressed in 
| J | units. 

After each field update an avalanche of spin flips propagates until the en- 
ergy is minimized for the applied field value. The distribution of avalanches 
is integrated along the hysteresis loop, and averaged over N run networks of 
given kind and size. This method has been applied both in simulations [22] 
and experiment in disordered magnetic systems [23]. The underlying physi- 
cal mechanism in these systems is the pinning of domain walls by quenched 
disorder of the network structure. It leads to the Barkhausen noise of flip- 
ping of magnetic domains [24]. The role of disorder in the fractal nature of 
the observed Barkhausen noise has not been fully understood [25]. The idea 
that the underlying processes of self-organized criticality [26] are responsible 
for the occurrence of the scale-free distributions of avalanches could not be 
proved in the general case. Often the observed finite avalanche cut-offs can be 
related to the finite sizes of the domains. Here we would like to stress that 
the controversy about the domain sizes does not appear in our networks. The 
only relevant scale is the network size. Likewise, the concept of domain walls 
has a different meaning, as discussed below. In our numerical experiment we 
can make distinction between the actual change of the magnetization A m and 
the number s of flipped spins within an avalanche s, which can be larger than 
A m (see Fig 5). 

The results for the distribution of avalanches N s (s) for different kinds of net- 
works are shown in Fig. 6. The plots are composed of two parts: For small 
values of the avalanche size s, the log-log plot is an approximately straight 
line; For larger sizes s we observe some oscillations with increasing amplitude. 
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Fig. 4. The field dependence of (a) the total number of flipped spins s(H ) /N in given 
field H and (b) density of walls (the number of antiferromagnetic pairs N\pp(H)) 
normalized to the total number of bonds Np for various kind of networks. The field 
is expressed in \J\ units. 

In the case of trees (M = 1, absence of frustration) the plots end with a well- 
defined maximum, corresponding to the largest avalanche along the hysteresis 
loop. The size of this avalanche is larger in the scale-free tree, compared with 
the exponential tree of the same number of links. For larger values of M, the 
multipeaky character of the plots N s (s) vanishes. However, for small s we ob- 
serve the power-law behavior of the spectra N s (s) in the whole investigated 
range of M, i.e. from 1 to 25. The values of the effective exponent a — defined 
by the relation N s (s) oc s~~ a — is constant for the exponential networks and 
it slightly increases with M for the scale-free networks [32]. 
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Fig. 5. Ratio of the magnetization change A m and the number s of the flipped 
spins within an avalanche. 



The distribution of the time r of duration of avalanches is shown in Fig. 7. 
Here, one time step is denned as the time necessary to scan all N spins in 
the system — spin by spin — to check if they flip or not. The time of the 
avalanche duration r is defined as the number of time steps which is necessary 
to get a stable state, i.e., when all flips die out at a given field H. As we 
see, the range of avalanche durations r covers only one order of magnitude. 
The reason is the small- world effect [3], i.e., short distance between the nodes 
in random graphs. As a consequence, the obtained range of r increases with 
the network size N only as log(iV). For the networks investigated here the 
avalanches never take more than nine or ten time steps. The diameters i [33] 
of the networks considered in this work, shown in Tab. 1, are comparable to 
the average avalanche duration. 

The data in Fig. 4(a) for SFN are compared with those in smaller graphs for 
M — 1 [Fig. 8(a)] and M = 2 [Fig. 8(b)]. In these figures the field is normalized 
to .f/max, which is different for different graph size, i.e., -ff max ~ 100, 100,200 
and 300 for N = 10 3 , 2 ■ 10 3 , 5 ■ 10 3 and 10 4 , respectively. The logarithmic plots 
in Fig. 8 reveal that the maximal avalanches normalized by number of nodes 
iV do not decay with the system size N. This means that the largest jump 
in the magnetization on the hysteresis curve persists when the network size 
increases, suggesting that true criticality of the hysteresis loop [27] occurs in 
the scale-free networks. Note that for trees, there is a gap of the avalanche 
spectrum between H = i? max and H = 0. This gap vanishes for M = 2. This 
is visible also in Fig. 6, where the largest avalanches form sharp and separate 
maxima of the curve N s (s). 
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(a) scale-free, M=1, N=10 4 , N^IO 3 



(b) exponential, M-1, N=10 4 , N run =10 3 
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(c) scale-free, M-5, N=10 4 , N run =10 3 
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(e) scale-free, M.15, N=10 4 , N run =10 3 
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(d) exponential, M-5, N=10 4 , N run =10 3 
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Fig. 6. Histograms of size s of avalanches during hysteresis experiment. 

When temperature increases, the hysteresis loops observed in the numerical 
experiments become more noisy, their area decreases and finally vanishes for 
T ~ 0.5 [| J\/k B \- Two examples of such loops obtained with Metropolis al- 
gorithms [28] are shown in Fig. 9 for T = 0.05 and T = 0.1 [\J\/k B ] [29]. 
The avalanche size distribution N s (s) preserves its power-law character for 
the above values of temperature. However, the thermal noise does not allow 
to separate sharp peaks for large s, which are visible in the case of T = (Fig. 
6). We note, that the calculations for T > are time-consuming. Still, we can 
deduce from the results of the simulation that the magnetization curves of the 
spin- networks at finite temperature are similar to those at T = 0, as long as 
the thermal energy is much smaller that the energy of the node-node magnetic 
interaction between nearest neighbors. 

A typical running for iV = 10 4 and M = 2 takes about half of minute of 
CPU time on Intel Itanium 2 (1.5 GHz) processor. For statistics given here 
(Arun = 10 3 ) it gives over eight hours of CPU time per one network and one 
branch of hysteresis m(H). 
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Fig. 7. Histograms of duration time r during hysteresis experiment. 
5 Discussion and conclusions 



The results show that the connectivity distribution in complex networks with 
antiferromagnetic bonds leads to broad distributions of avalanches. In the 
scale-free networks, the maximal degree is larger than in the exponential net- 
works of the same size. This difference produces broader loops for SFN than 
for EXN. On the other hand, for networks with larger connectivity (M = 2) 
the right parts of histograms of avalanches are irregular, whereas for trees 
(M = 1) we observe fine structure of sharp maxima with zero signal between 
them. 

It is possible to interpret the size of maximal avalanches s max in terms of the 
structure of the investigated networks. It appears that s max is comparable with 
the number of nodes with minimal degree. On the other hand, s max happens 
to be equal to the number of spins which flip at field between H = —M + 5 
and H = — M — 1 + 8 for decreasing field H. Let us trace these rules in detail. 

In Fig. 6 we see that maximal avalanches for SFN trees contain approximately 
7100 flipping spins, i.e. 71 percent of the system size. This value corresponds 
to the height 0.72 of the maximum in Fig. 4(a), which happens to appear at 
H = —\J\ again for the decreasing field branch. On the other hand, the degree 
distribution for SFN is known to be [30] 



P(k > M) 



2M(M + 1) 



(2) 



fc(Jfe + l)(Jfe + 2) 

For trees, the minimal degree is k min = M = 1, then, for SFN we get 
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(a) scale-free tree (M=1): N run =1 
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(b) scale-free graph (M=2): N run =1 
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Fig. 8. Fraction of flipped spins s/N normalized to the network size N in increasing 
field H for various values N in the scale-free networks with (a) M = 1 and (b) 
M = 2. 



P(k = 1) = 2/3. These three numbers, in first column of numbers in Ta- 
ble 1 (scale-free, M=l) are close to each other; our interpretation is that they 
give approximate information on the same process. 

The same expression (2) applies to SFN with M — 2, i.e. k min = 2, as the 
degree cannot be smaller than M. Calculating P(k = 2) from Eq. (2) gives 
P(k = 2) = 1/2. This is to be compared with the size of the maximal avalanche 
for SFN, M = 2 in Fig. 6, and to the fraction of spins which flip at field close 
to H = —2\J\, as shown in Fig. 4(a). The numerical values are given in the 
second column (scale-free, M — 2) of Tab. 1. As we see, they all of them are 
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Fig. 9. (Color on-line). Hysteresis loop for finite temperature T = 0.05 and 0.1 
[|J|//cb]. The field H is measured in |J|. 

Table 1 

Fraction of spins flipped at H = —M [Fig. 4(a)] and maximal avalanche size s max /iV 
(Fig. 6) as compared with fraction of spins with degree equal to M [4,17,30] for 
various networks. The effective exponents a and (3 of power law distributions are 
included for N s (s) and N m (A m ), respectively. The last line shows network diameters 
I [33] for N = 10 4 . 

scale-free exponential 



M 


i 


2 


1 


2 


s for H = -M 


0.72 


0.49 


0.41 


0.23 




0.73 


0.54 


0.44 


0.25 


P(k = M) 


2/3 


1/2 


1/2 


1/3 


a 


1.38 


1.39 


0.97 


1.00 


P 


1.56 


1.48 


0.92 


0.93 


e 


9.1 


5.1 


15.6 


5.1 



close to 0.5. Again, this means that the largest avalanche contains spins at all 
nodes with M — 2, and this avalanche happens also at H = —\J\M. 

The same evaluation for EXN, i.e., 0.41 for M — 1, and 0.23 for M = 2 from 
Fig. 4(a), when compared with N s (s) gives a qualitative accordance. As a 
rule, the size of maximal avalanches is comparable with the system size. This 
results suggests that the self-organized criticality (SOC) might be a possible 
mechanism of the spin-reversal dynamics in the antiferromagnetic disorder- 
free networks. We note that for EXN trees, the number of leaves of minimal 
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(a) scale-free, M=1,N=1CP, N ru 



(b) exponential, M.2, N=10 4 , N run =10 3 
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Fig. 10. Histograms of changes of magnetization A m of avalanches during hysteresis 
experiment. 

degree k = 1 can be evaluated from P(k > M) — 2~ k [4] as 1/2 for M = 1. For 
M = 2, the distribution P(k > M) — 3/4 • (3/2) ~ k [17] gives P(k = 2) = 1/3. 
Comparison of the results for EXN is included in third and fourth column 
of Tab. 1. These numbers prove, that in the hysteresis experiment, maximal 
avalanches for the descending field appear when H = — k niin \ J\ = —M\J\. The 
spins flipping in these avalanches occupy the nodes with minimal degree M. 
The observed quantitative differences might be attributed to the finite net- 
work sizes in the numerical calculations; this produces deviations of simulation 
results from theoretical formula, which are exact in the thermodynamic limit. 



Left parts of the histograms, i.e. the spectra for small avalanche sizes, are 
straight lines in the log-log scale. This means that N s (s) oc s~ a , with a close 
to 1.00 for the exponential networks and a slightly increasing with M for 
the scale-free networks. It is not clear if this variation of a for the scale-free 
networks is a consequence of a M-dependent finite size effect, or it means 
just a lack of universality. We cannot verify it, because for large networks, 
the numerical algorithm becomes very time-consuming. However, this scaling 
behavior is valid only in the range of the avalanche size s not much larger 
than one order of magnitude [11]. The latter apply also to the distribution 
N(A m ) oc A m @ (see Fig. 10). The effective exponents a and (3 are also shown 
in Tab. 1. Spins which flip in this way are neither leaves (for M = 1), nor 
"leaves with two stems" (for M — 2), but they form "central" (here the 
oldest) parts of the networks. The results suggest that the spin flips in the 
better connected parts of the network form small avalanches with power-law 
distributions. 



Basing on these results we expect that avalanches and critical behavior may 
occur in various systems with broad connectivity distributions where some 
negative feedback exists between local degrees of freedom. In our spin-networks 
it is the antiferromagnetic character of interaction which causes this kind of 
feedback. More direct interpretation depends on particular systems. In par- 
ticular, avalanches seem possible in complex magnetic systems as quantum 
magnetic cellular automata with magnetostatic interaction between particu- 
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lar monodomain elements [31]. The size of avalanches is limited only by the 
size of the network. Then, this kind of collective flipping is related to an insta- 
bility of the system. Design of an appropriate spatial structure, in particular 
by increasing clustering of the networks [11] allows to tune and control the 
magnetization curve. 

Concluding, our numerical results indicate that the positions of the largest 
peaks in the spectra of avalanches, which measure the size of the largest 
avalanches, agree well with the numbers of spins flipped at the field equal to 
the connectivity parameter M. The latter can be interpreted as the numbers 
of nodes with minimal degree for networks with different topologies and con- 
nectivities. These results agree qualitatively with known theoretical formula 
for the degree distributions. 
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